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HE: 全 球 卫 星 导 航 系统 多 系统 选 星 问题 能 够 转化 为 有 约束 条 件 的 多 目标 优化 问题 求解 ， 可 以 同时 优化 几何 精度 因子 
(GDOP) 和 选 星 数目 这 两 个 目标 ， 从 而 在 减少 接收 机 运算 量 的 同时 获得 良好 的 定位 精度 。 本 文 提 出 了 种 基于 NSGA-I 算 
法 的 多 目标 快速 选 星 方法 ， 该 方法 利用 选 星 问题 的 序 贯 性 生成 初始 种 群 ， 改 进 了 约束 处 理 方法 ， 并 选取 合 oo 
函数 做 出 选 星 决策 。 通 过 仿真 ， 证 实 该 选 星 方法 具有 良好 的 可 靠 性 和 实时 性 ， 且 不 依赖 于 卫星 的 几何 位 置 分 布 ， T 
障碍 或 者 遮挡 的 复杂 情况 。 
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全 球 卫 星 导 航 系统 (Global Navigation Satellite System, GNSS) 发 展 迅 猛 ， 从 前 美国 全 球 定 位 系统 一 
家 独 大 的 局 面 已 被 打破 ， 随 着 俄罗斯 GLONASS 系统 、 欧 洲 伽 俐 略 〈Galileo) 系统 和 中 国 北斗 (BDS) 系 
统 的 先后 崛起 ， 四 大 卫星 导航 系统 并 立 共 存 的 格局 已 基本 形成 。 多 系统 组 合 导航 已 成 为 今后 的 主要 发 展 趋 
势 ， 为 用 户 提供 更 为 精确 可 靠 的 定位 、 导航 和 授时 服务 na。 3。 现 今 ， 即 使 是 低 成 本 的 接收 机 都 可 以 同时 接收 
r 多 个 全 球 卫星 导航 系统 的 信号 ， 但 是 过 多 的 元 余 信 息 会 在 增 大 定位 解 算 的 计算 量 的 同时 增 大 功 耗 ， 严 重 影 
a Mel ECA LAE PEA SEIN E, E, ae BE A OT PAF a WL EMR RE ATT ERE, BAERE MEN 
LO 选 星 是 选取 最 优 的 可 见 卫 星 组 合 的 过 程 ， 通 常 利 用 几何 精度 因子 C Geometric Dilution of 
aa Precision, GDOP) NARRE NFU NE, FCP ASCE E ce ae A i IE. Be 
© 导航 系统 选 星 算法 是 在 原 有 的 单 系统 选 星 算 法 的 基础 上 发 展 起 来 的 ， 例 如 基于 卫星 几何 构 型 的 最 大 凸 多 面 
本 体积 法 外 和 基于 高 度 角 和 方位 角 的 选 星 算法 %9 等 。 文 [7] 提 出 了 一 种 适用 于 典型 的 城市 障碍 环境 的 选 星 算 
法 ， 通 过 单 系统 选 星 构成 初始 集合 ， 然 后 逐步 扩大 选 星 范围 来 获得 有 最 小 几何 精度 因子 值 的 选 星 结果 。 文 
[8] 通 过 计算 每 颗 卫 星 对 几何 精度 因子 的 贡献 ， 然 后 剔除 贡献 最 小 的 卫星 达到 选 星 目的 。 进 化 算法 和 机 器 学 
习 也 被 大 量 的 运用 到 了 选 星 过 程 中 ， 包 括 粒 子 群 算法 外、 模拟 退火 算法 0、 神 经 网 络 算法 50、 遗 传 算法 呈 - 
41 SHARPLES. KH, 文 [9] 认 为 由 于 神经 网 络 算法 需要 大 量 的 训练 过 程 ， 进化 算法 在 处 理 选 星 
问题 时 的 速度 更 快 。 文 [14] 通 过 比较 一 系列 选 星 方法 ， 结 果 表 明基 于 支持 向 量 机 和 遗传 算法 的 选 星 性 能 
s 优 。 
ee 上 述 的 选 星 方法 大 多 是 在 先 验 性 的 确定 选 星 数目 的 条 件 下 ， 选 择 几 何 精度 因子 最 优 的 可 见 卫 星 组 合 ， 
‘S 因此 都 是 单 目 标 优化 方法 。 而 在 运用 多 系统 导航 的 今天 ， 任 意 地 点 的 平均 可 见 导航 卫星 数目 将 超过 40 颗 ， 


= 选 多 少 颗 星 用 于 接收 机 的 定位 解 算 ， 也 成 为 了 愈加 重要 的 考量 因素 。 在 这 方面 ， 以 往 的 选 星 算法 不 能 依据 
"一 卫星 几何 布局 的 优 劣 ， 机 动 地 确定 选 星 数目 ， 有 着 固有 的 局 限 性 。 当 选 星 数目 过 少时 ， 几 何 精 度 因 子 值 增 
-= 加 从 而 影响 定位 精度 ， 而 选 星 数目 过 多 时 ， 几 何 精度 因子 值 减 小 但 是 大 大 增加 接收 机 的 计算 负荷 和 功 耗 ， 
几何 精度 因子 和 选 星 数目 是 两 个 改善 方向 相互 制约 的 目标 ， 其 中 一 个 目标 的 优化 必然 致使 另 一 个 目标 的 性 
能 降低 。 因 此 ， 选 星 问题 可 以 转化 为 多 目标 优化 问题 求解 。 这 样 可 以 减少 以 往 选 星 算法 中 ， 先 验 生 的 确定 
选 星 数目 所 带 来 的 主观 影响 ， 同 时 不 再 拘泥 于 固定 选 星 数 带 来 的 对 几何 精度 因子 极 值 的 约束 ， 提 高 了 选 星 
的 机 动 性 。 这 种 多 目标 选 星 方式 ， 能 够 得 到 几何 精度 因子 和 选 星 数目 这 两 个 目标 进行 权衡 之 后 的 最 优选 星 
方案 ， 可 以 更 好 地 满足 接收 机 定位 解 算 的 精 度 和 实时 性 需求 。 
本 文 开 创 性 地 将 多 卫星 导航 系统 选 星 问题 作为 有 约束 条 件 的 多 目标 优化 问题 进行 研究 ， 并 给 出 了 数学 
描述 。 进 一 步 提出 了 一 种 基于 快速 非 支 配 排序 遗传 算法 (Non-dominated Sorting Genetic Algorithm- 
I，NSGA-I) 的 快速 选 星 方法 ， 能 在 同时 优化 几何 精度 因子 和 选 星 数目 这 两 个 目标 ， 并 在 各 种 全 球 卫 星 导 
航 系 统 组 合 条 件 下 进行 仿真 和 分 析 。 结 果 验 证 了 该 方法 的 解决 多 系统 选 星 问题 的 有 效 性 和 实用 性 。 


1 多 目标 选 星 问题 的 数学 描述 
多 目标 选 星 问 题 是 选择 最 优 的 卫星 组 合 的 过 程 ， 目 的 是 同时 得 到 最 小 的 选 星 数 目 和 几何 精度 因子 数值 。 
为 了 能 将 选 星 问题 转化 为 数学 问题 ， 以 便 采用 合适 的 算法 求解 ， 同 时 降低 选 星 求解 的 运算 量 ， 本 文 依托 先 
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星 问 题 的 实际 特性 ， 将 数学 模型 表述 为 
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. (1) 
RTA] BX = ( xixzxs … Xn)" 代表 不 同 的 选 星 方案 ， 所 在 的 空间 为 决策 空间 ，n 为 可 见 导航 卫星 的 总 数 。 
Rh, AREJA LP, x = 1 表示 卫星 被 先入 ，%y = 0 表示 卫星 未 入 选 ， 这 样 的 设计 可 以 方便 地 对 选 
星 方案 进行 二 进 制 编码 。 两 个 目标 函数 fx) 和 fx) 所 在 的 空间 为 目标 空间 ， 分 别 代 表 选 星 方案 和 的 几何 精 
度 因子 值 和 选 星 数目 。gi(x) 表 示 约 束 条 件 ， 确 定 了 解 的 可 行 域 。 对 于 个 卫星 导航 系统 而 言 ， 定 位 求解 时 
需要 K+3 颗 可 见 卫 星 ， 监测 还 有 故障 检测 和 排除 的 要 求 ， 则 需要 选择 更 多 的 卫星 。 
通常 情况 下 ， 定 位 时 选用 的 卫星 越 多 ， 定 位 精度 越 高 ， 然 而 ， 考 虑 到 接收 机 的 解 算 能 力 ， 需 要 设置 选 星 数 
目的 上 限 mmo PAE, AE BRDU 件 的 多 目标 优化 问题 求解 。 
多 系统 定位 需要 考虑 到 系统 间 的 不 一 致 性 。 在 目前 所 能 达到 的 测量 精度 水 平 下 ， 系 统 间 坐标 参考 框架 
的 偏差 可 以 忽略 不 计 ， 所 以 需要 重点 考虑 系统 间 的 时 间 偏 差 。 对 于 用 户 而 言 ， 时 间 偏 差 可 以 并 入 接收 机 钟 
差 ， 每 增加 一 个 系统 ， 就 增加 一 个 钟 差 参 数 ， 观 测 和 矩阵 态 也 需要 进行 相应 的 扩展 。 因 此 ， 多 系统 的 几何 精 
度 因子 和 测量 矩阵 鼠 可 以 表示 为 5 
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(2) 
(3) 
y 其 中 ， 下 标 y=, B, C … 表 示 不 同 的 导航 系统 ， 例 如 GPS. GLONASS, Galileo 和 | 
2 PEM n=nytngtinct....0 在 观测 矩阵 及 中 ， 孔 是 一 个 力 x3 维 的 方向 余弦 和 矩阵 ， 分 别 为 力 维 列 
wW ”向 量 ， 用 于 解 算 系 统 钟 差 。 
2 快速 选 星 方法 


多 目标 优化 旨 在 获得 一 组 多 目标 函数 折衷 解 的 集合 ， 称 为 Pareto 最 优 解 集 或 者 非 支 配 解 集 ， 然 后 通过 
一 定 的 决策 规则 找到 其 中 最 满意 的 方案 0。 遗 传 算法 (Genetic Algorithms, GA) 作为 一 种 模拟 自然 选择 的 
启发 式 随 机 搜索 算法 ， 被 广泛 运用 于 多 目标 优化 领域 。 在 现 有 的 多 目标 遗传 算法 中 ，NSGA-I 算 法 求解 多 
- 系统 选 星 问题 最 为 合适 ， 该 算法 运用 了 精英 保留 机 制 、 快 速 非 支配 排序 方法 和 计算 拥挤 距离 的 多 样 4 性 保持 
© 策略 3。NSGA-II 算 法 的 基本 思想 为 首先 产生 一 个 规模 为 N 的 初始 种 群 ， 然 后 在 每 一 代 的 进化 中 ， 父 代 种 
N BE PAI AR AN EER, LMR 3 个 基本 操作 ， 得 到 子 代 种 群 0,， 将 父 代 种 群 和 子 代 种 群 合 并 P,U 0, 之 
> 后 ， 进 行 快速 非 支配 排序 ， 再 通过 拥挤 度 计算 确定 个 体 虚拟 的 适应 度 ; 最 后 选择 最 好 的 NN 个 个 体 生成 新 的 
PEE Pris 依次 循环 ， 直 到 满足 结束 条 件 。 
由 此 可 知 ，NSGA-II 算法 的 主要 优点 是 快速 非 支配 排序 方法 能 降低 计算 的 复杂 度 ， 并 且 不 用 外 部 种 群 
储存 非 支 配 个 体 ， 因 此 初始 种 群 的 规模 极 大 地 影响 算法 的 性 能 ， 需 要 进行 重点 分 析 ?3。 同 时 ，NSGA-II 算 
法 中 使 用 的 遗传 算 子 ， 包 括 模拟 二 进 制 交叉 算 子 〈Simulated Binary Crossover) 和 多 项 式 变异 算 子 
(Polynomial Mutation) ， 并 不 适用 于 二 进 制 编码 的 选 星 问题 ， 因 而 需要 设计 更 为 合适 的 遗传 算 子 。 此 外 ， 
获得 Pareto 最 优 解 集 只 是 解决 多 目标 优化 问题 的 第 1 步 ， 之 后 同样 重要 的 步骤 是 通过 一 定 的 决策 规则 从 
Pareto 最 优 解 集中 选择 一 个 满意 解 。 因 此 ， 依 据 NSGA-II 算法 的 上 述 特 性 ， 本 文 提 出 了 一 种 改进 的 基于 
NSGA-I 算法 的 多 目标 选 星 方法 ， 整 体 思 路 沿用 NSGA-I 算 法， 但 是 在 生成 初始 种 群 时 ， 利 用 了 选 星 问题 
的 序 贯 性 ， 改 进 了 处 理 约束 方法 ， 同 时 选取 合适 的 遗传 算 子 ， 并 设计 最 后 决策 时 使 用 的 效用 函数 。 


2. 1 初始 化 种 群 


多 目标 选 星 问题 可 以 方便 地 采用 二 进 制 编码 方式 。 假 设 某 时 刻 的 可 见 卫 星 总 数 为 na， 按照 顺序 编排 之 
后 ， 将 每 颗 卫 星 都 设 为 一 个 基因 ， 入 选 时 记 为 1， 未 入 选 时 记 为 0。 编 码 后 的 所 有 可 见 卫 星 组 成 一 个 长 度 为 
n 的 二 进 制 字符 串 ， 称 为 一 条 染色 体 。 每 条 染色 体 代 表 一 个 可 能 的 选 星 方案 。 通 常情 况 下 ， 初 始 种 群 由 一 
组 随机 生成 的 染色 体 组 成 。 
连续 选 星 问 题 可 以 看 作 是 一 个 序 贯 决策 (Sequential Decision-Making) 过 程 ， 需 要 在 动态 环境 下 ， 序 
贯 性 地 在 每 个 时 刻 ， 根 据 所 处 的 状态 和 以 前 的 记录 ， 从 一 组 可 行 方案 中 做 出 最 优 决 策 。 文 [19] 进 过 详尽 的 
分 析 ， 证 实 了 每 一 个 时 刻 的 最 优选 星 组 合同 样 也 是 一 定时 间 间 隔 范 围 内 的 最 优选 星 组 合 ， 当 间隔 时 间 设 置 
为 1min 时 ，95% 的 时 间 内 可 以 满足 这 个 规律 。 图 1 对 比 展 示 了 一 个 全 球 定位 系统 和 北斗 导航 系统 组 合 在 
相 邻 的 前 后 两 个 时 刻 的 最 优选 星 方案 ， 仰 角 门 限 为 5?， 选 星 数 目 为 7。 由 图 可 知 ， 在 这 个 时 间 跨 度 内 ， 最 
优选 星 方案 中 有 5 颗 可 见 卫 星 依然 匹配 ， 有 2 颗 发 生 了 变化 。 发 生变 化 是 因为 随 着 时 间 的 推移 ， 有 卫星 出 


现 或 者 消失 在 视野 范围 内 ， 这 种 情况 也 仅仅 发 生 在 5% 的 时 间 内 。 结 果 证 明 ， 可 以 利用 上 一 时 刻 的 选 星 结 
果 生 成 下 一 时 刻 的 选 星 方案 。 


图 1 全 球 定位 系统 和 北斗 导航 系统 前 后 两 个 时 刻 的 最 优选 星 方案 对 比 
Fig. 1Comparison of the GPS/BDS optimal satellite subsetsat two time points 


基于 选 星 问 题 的 序 贯 性 ， 前 一 时 刻 的 Pareto 最 优 解 集 和 下 一 时 刻 的 最 优 解 集 ， 应 该 相同 或 者 具有 极 大 的 
相似 性 。 因 此 可 以 将 前 一 时 刻 的 Pareto 最 优 解 集 译 码 后 ， 获 得 最 优 的 选 星 结果 ， 然 后 去 掉 其 中 下 一 时 刻 不 
可 见 的 卫星 ， 生 成 新 的 种 群 ， 作 为 下 一 时 刻 的 初始 种 群 P。。 这 样 就 可 以 利用 前 一 刻 的 Pareto 最 优 解 集 ， 生 
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成 下 一 时 刻 的 初始 种 群 ， 具 体 过 程 如 图 2。 
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图 2 二 进 制 编码 和 基于 序 贯 性 生成 初始 种 群 


Fig.2Binary encoding and Population initialization with sequential decision method 


基于 上 述 设置 ， 初 始 种 群 已 经 接近 最 优 解 集 ， 于 是 利用 较 小 的 种 群 规 模 ， 经 过 较 少 的 进化 代数 ， 就 可 
以 得 到 满意 的 选 星 结果 。 这 将 大 大 减少 计算 的 复杂 度 ， 并 且 能 提高 连续 选 星 的 效率 ， 这 样 的 改进 方式 也 可 
以 广泛 运用 于 其 他 的 选 星 方法 中 。 


2.2 约束 处 理 
有 约束 条 件 的 多 目标 优化 算法 的 收敛 性 和 分 布 性 的 优 和 劣 ， 取 决 于 约束 处 理 技术 和 多 样 性 维护 策略 ， 现 
有 的 约束 处 理 技术 包括 惩罚 函数 法 、 随 机 排序 法 、 可 行 性 准则 法 和 混合 算法 等 外。 其 中 惩罚 函数 法 需要 设 


置 合适 惩罚 系数 ， 不 然 就 会 早熟 或 者 停滞 ， 随 机 排序 法 稳定 性 较 差 ， 而 基于 解 之 间 的 约束 支配 关系 ， 可 行 


性 准则 认为 后 : (可 行 解 支配 不 可 行 解 ，(2) 可 行 解 之 间 适 应 度 高 的 解 支 配 适应 度 低 的 解 ，(3) 不 可 行 解 之 


间 约 束 违反 度 小 的 解 支配 约束 违反 度 大 的 解 。 该 算法 由 于 过 分 强调 了 可 行 解 优 于 不 可 行 解 ， 容 易 陷 入 局 部 


最 优 。 


可 行 解 所 占 比例 的 阔 值 ， 是 个 单调 递减 函数 ， 其 表达 式 为 ; 


a=-0.35，p=0.2。 在 设 定 上 述 阔 值 的 基础 上 ， 生 成 新 种 群 的 Pi 的 选择 操作 如 下 : 


而 选 星 问题 的 最 优 解 往往 存在 于 约束 边界 ， 为 了 保持 解 的 多 样 性 ， 约 束 违反 度 较 小 的 不 可 行 解 也 需要 
在 种 群 中 得 到 一 定 的 保留 。 为 此 ， 本 文 提 出 了 一 个 保持 一 定 比 例 不 可 行 解 的 新 的 选择 算 子 p, 代表 种 群 中 不 


ER, t 为 进化 代数 ;Giwwx 为 最 大 进化 代数 ，a 和 45 均 为 [0,1] 范 围 内 的 可 调 参 数 ， 在 选 星 问题 中 可 以 设置 为 


(1) 合并 规模 为 NN 的 父 代 和 子 代 种 群 ， 生 成 规模 为 2N 的 REPU Os 设置 新 种 群 Po= 纪 i=1. 


( 2 ) 如 果 R, 中 不 可 行 解 的 个 数 Ninfeasible < int(pN) ? 选择 所 有 的 不 可 行 解 Xinfeasible 5 如 R Ninteasibezint(oN) ? 
则 利用 公式 ， 求 不 可 行 解 的 约束 违反 程度 ， 选 择 约束 违反 程度 较 小 的 intloM) 个 不 可 行 解 xinteasinie; Pei= 


Put U Xinfeasible o 


(3) 对 玉 ,中 的 可 行 解 进行 快速 非 支 配 排序 ， 并 标记 非 支配 等 级 Ff， 二 1,2,3... 当 |Piml+|FI<N 时， 按照 非 
支配 等 级 高 低 依次 选择 个 体 PaPa UF, i=i+1; 然后 依据 拥挤 度 排序 ， 在 非 支配 等 级 Ff 中 选择 分 布 性 较 


好 的 可 行 解 填充 ， P= Pr UF[l: N-Pnil]。 


改进 的 约束 处 理 方式 通过 保留 一 定数 量 的 不 可 行 解 ， 增 强 了 在 可 行 域 边缘 的 搜索 能 力 ， 一 定 程度 上 维 


护 了 种 群 的 多 样 性 ， 避 免 了 早熟 现象 的 发 生 。 
2. 3 遗传 算 子 


遗传 算 子 包括 选择 算 子 、 交 叉 算 子 和 变异 算 子 。 变 异 操作 指 以 一 定 的 概率 改变 个 体 编 码 串 中 的 某 个 或 
者 某 些 基因 位 置 的 过 程 。 对 于 选 星 问题 可 以 采用 基本 位 变 腊 方法 ， 依 据 变异 概率 指定 变异 点 ， 然 后 对 指定 


的 变异 点 上 的 基因 值 做 取 反 运算 即 可 。 交 叉 算 子 的 选用 一 般 与 编码 方式 有 关 ， 


主要 有 单 点 交叉 、 多 点 交叉 、 均 匀 交 叉 、 循 环 交 叉 和 部 分 匹配 交叉 等 。 为 了 


用 异 或 交叉 算 子 ， 采 用 轮 盘 法 则 ， 随 机 选择 父 代 进行 交叉 操作 ， 其 操作 方式 如 
fi folififofifo}i| 
ofotifofofifi}y| 


1} 11 
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OT oT 


图 3 异 或 交叉 算 子 
Fig.3XOR crossover 


2.4 效用 函数 


二 位 进 编码 常用 的 交叉 操作 


效 保留 优秀 的 基因 ， 本 文选 


图 3。 


多 目标 进化 算法 通常 分 为 优化 和 决策 两 个 步 又 ， 所 采用 的 技术 一 般 分 为 3 类 : 先 验 决策 方法 (决策 一 


优化 )、 后 验 决 策 方法 (优化 一 决策 ) 和 交互 决策 方法 〈 优 化 一 一 决策 ) I, 


本 文采 用 后 验 决 策 方法 ， 在 


利用 改进 的 NSGA-I 算 法 获得 选 星 问题 的 Pareto 最 优 解 集 之 后 ， 用 一 个 更 高 层次 的 效用 函数 选择 一 个 偏好 


解 。 针 对 双 目 标 优化 的 选 星 问 题 ， 设 计 的 效用 函数 表示 为 阶段 
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即 为 选 星 结果 。 
这 种 后 验 决 策 方法 ， 和 有 
特性 ， 并 行 处 理 各 个 目标 ， 


结 


根据 上 述 


其 中 ， 为 Pareto 最 优 解 集中 的 第 7 个 选 星 方案 ， 两 个 目标 函数 的 最 大 值 和 最 小 值 分 
的 权 值 。 利 用 效用 函数 计算 得 到 Pareto 最 优 解 集中 每 个 选 星 方 
经 过 解码 后 ， 可 利用 选 星 结 


果 进 行 定位 解 算 。 
E 够 有 效 地 控制 决策 者 的 喜好 带 来 的 主观 影响 ， 
在 具备 全 局 最 优 搜索 性 的 条 件 下 ， 
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别 为 和 » Wi 和 W2 为 设计 
适应 值 最 小 的 选 星 方案 最 优 ， 


案 的 适应 值 后 ， 


并 且 能 够 很 好 地 利用 遗传 算法 的 
效 地 优 中 选 优 ， 最 终 得 到 满意 的 选 星 


快速 高 


术 分 析 和 改进 方案 ， 基 于 NSGA-II 算法 的 多 目标 快速 选 星 方法 的 流程 图 如 图 4。 


yi 
初始 化 种 群 


遗传 操作 选择 、 


交叉 和 


变异 ) 后 ， 产 生子 代 


子 代 和 父 代 合并 


计算 种 群 中 各 个 选 星 
方案 的 目标 函数 值 和 


约束 违反 程度 


非 支配 排序 和 拥挤 度 计 


算 ， 选 择 产 生 新 种 群 


利用 效用 函数 决策 
输出 选 星 结果 


图 4 基于 NSGA-I 算法 的 多 目标 快速 选 星 方法 的 流程 图 
Fig. 4 Flow chart of fast satellite selection method based on NSGA-II 


3 仿真 结果 和 实验 分 析 


为 了 测试 基于 NSGA-I 算 法 的 多 目标 快速 选 星 方法 的 性 
导航 系统 全 面 运行 时 的 情况 


境 为 各 个 卫星 

据 实际 情况 ， 采 用 了 真实 的 星 历 。 对 于 中 
星 、3 颗 倾 斜 地 球 同步 轨道 卫星 

GPS+BDS 和 GPS+BDS+GLO。 


能 ， 本 文 用 自 编 软件 进行 了 仿真 分 析 。 念 真 环 


。 其 中 ， 全 球 定位 系统 和 GLONASS 由 于 已 经 部 署 完成 ， 于 是 根 


ESB 


中 包括 5 颗 地 球 同步 轨道 


上 斗 系统 按照 35 MELT HA, E 


和 25 颗 中 地 球 轨道 卫星 。 依 据 普 遍 性 ， 采 用 了 两 种 组 合 方案 : 
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Ca) HYD PB Cb) 几何 精度 因子 
图 5 GPS+BDS 选 星 前 可 见 卫星 数目 和 几何 精度 因子 
Fig. 5 (a) The number of visible satellites of GPS+BDS and (b) GDOP prior to the selection process 
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(a) 选 星 数目 Cb) 几何 精度 因子 
图 6GPS+BDS 选 星 数目 和 选 星 后 的 几何 精度 因子 
Fig. 6 (a) The number of selected satellites of GPS+BDS and (b) GDOP following the selection process 


5 AANA GPS+BDS 组 合 系统 在 选 星 前 ， 全 球 范围 内 的 可 见 卫 星 数 和 几何 精度 因子 ， 分 辨 率 为 
5?x5"， 高 度 截 止 角 设 为 5"。 图 6 是 该 时 刻 的 选 星 数目 和 选 星 后 的 几何 精度 因子 ， 设 选 星 数目 的 最 大 值 
mmax=15。 选 星 前 后 的 具体 结果 见 表 1。 由 此 可 知 ， 在 该 时 刻 全 球 的 平均 可 见 卫星 数 约 为 23 颗 ， 如 果 从 中 选 
取 7 颗 ， 将 有 =245157 种 组 合 方式 ， 而 可 用 的 卫星 数目 更 多 时 ， 选 星 的 组 合 方式 将 成 几何 级 数 增长 ， 这 对 
于 遍历 法 将 是 沉重 的 计算 负荷 ， 几 乎 不 能 实现 。 而 经 过 了 本 文 提出 的 选 星 过 程 ， 选 择 的 卫星 数目 均值 为 
10.5953， 为 选 星 前 的 46.7%， 使 用 的 卫星 数目 减少 ， 定 位 解 算 时 运算 量 的 可 以 减少 60%， 接 收 解 算 的 效率 
将 能 得 到 显著 提高 。 

在 选 星 前 的 几何 精度 因子 的 最 小 值 为 0.9652， 当 几何 精度 因子 大 于 1 时， 是 定位 误差 放大 因子 ， 当 几 
何 精度 因子 小 于 1 时 ， 说 明 卫 星 的 几何 布局 较 好 ， 于 是 具有 了 定位 误差 缩小 的 功能 。 在 选 星 后 的 几何 精度 
因子 最 大 值 不 超过 3， 依 然 可 用 性 良好 ， 而 且 选 星 后 的 几何 精度 因子 均值 为 1.5682， 仅 比 选 星 前 的 均值 增 
大 了 0.2 左 右 ， 说 明 选 星 方法 能 够 很 好 地 保持 接收 机 的 定位 精度 。 


表 1GPS+BDS 全 球 范围 选 星 前 后 统计 分 析 
Table 1 Statistical analysis of the number of selected GPS+BDS satellites and GDOP at the global scale before 
and after the satellite selection process 


选 星 数目 几何 精度 因子 
15 


| wmm | 15 | 32 | 227027 0.9652 1.9503 1.3168 


| wmm | 6 | 1 | 105953 1.2411 2.1662 1.5682 


在 GPS+BDS+GLO 系统 组 合 导 航 的 情况 下 ， 分别 对 选 星 问题 进行 仿真 ， 使 用 Intel Core i7-5500U 
(2.4GHz/L3 4M) 处 理 器 。 可 以 代表 仿真 时 间 为 24h， 采 样 间 隔 为 60s， 用 户 位 于 己 点 
(39.9°N, 116.3°E, 0) , 选 星 数目 的 最 大 值 Ma 设 为 该 时 刻 可 见 卫星 数目 的 60%。 高 度 角 的 设置 主要 是 为 
了 控制 接收 机 接收 卫星 的 数目 ， 通 常情 况 下 ， 由 于 树木 或 者 高 楼 的 遮挡 ， 以 及 多 路 径 效应 影响 ， 高 度 角 小 
的 卫星 信号 受到 的 干扰 也 会 比较 大 ， 定 位 解 算 的 结果 不 准确 ， 所 以 接收 机 需要 尽量 不 使 用 高 度 角 小 的 卫星 。 
在 实际 测量 中 ， 需 要 根据 周围 环境 的 开阔 程度 调整 高 度 截止 角 的 设置 ， 环 境 越 开 阔 ， 高 度 截 止 角 可 以 设置 
的 越 小 。 所 以 在 本 文 的 仿真 实验 中 ， 高 度 截止 角 分 别 设置 为 5*、15°* 和 30°*， 结 果 见 表 2 及 图 7。 


a 


K 2GPS+BDS+GLO 选 星 前 后 统计 分 析 
Table 2Statistical analysis of the number of selected GPS+BDS+GLO satellites and GDOP at the global scale 
before and after the satellite selection process 


| RBA 几何 精度 因子 


rice 
rans] | [om ian [eo 
Pee 


星 后 fame e he sr 


| [oe par | 750 | zie07 [a 


当 高 度 截止 角 为 5* 时 ， 选 星 数目 均值 为 选 星 前 卫星 数目 的 40%， 选 星 后 几何 精度 因子 的 最 大 值 小 于 
1.7， 与 选 星 前 的 最 大 值 相当 ， 同 时 选 星 后 的 几何 精度 因子 的 平均 值 ， 相 比 于 选 星 前 只 增加 了 0.26。 选 星 后 
的 几何 精度 因子 在 仿真 时 段 内 的 与 选 星 前 差别 不 大 ， 保 持 了 良好 的 可 用 性 。 当 高 度 截止 角 为 13" 时 ， 有 
部 分 时 间 段 内 ， 选 星 后 的 几何 精度 因子 值 比 选 星 前 的 还 要 小 ， 这 说 明 在 定位 解 算 的 时 候 确实 有 许多 宛 余 信 
息 ， 需 要 经 过 选 星 过 程 剔 除 ， 同 时 说 明 本 文 设计 的 选 星 方法 的 性 能 优良 ， 能 有 效 地 寻找 到 合适 的 选 星 方案 。 
这 种 优越 性 在 高 度 截止 角 为 30° 时 ， 表现 更 为 突出 ， 在 30% 的 时 间 段 内 ， 选 星 后 的 几何 精度 因子 都 优 于 选 
星 前 ， 这 也 说 明了 选 星 的 必要 性 。 由 图 7 可知 选 星 耗 时 分 布 在 10ms 至 70ms 之 间 ， 均 值 约 为 30ms， 这 证 明 
了 本 算法 的 实时 性 。 

由 于 本 文 提 出 的 选 星 方法 为 全 局 搜索 算法 ， 不 依赖 于 卫星 之 间 的 几何 分 布 ， 所 以 没有 适用 条 件 的 限制 ， 
特别 适用 于 有 遮挡 或 者 障碍 的 复杂 情况 。 同 时 ， 可 以 调整 算法 中 的 主要 参数 和 关键 步骤 ， 例 如 种 群 规 模 、 
进化 代数 、 约 束 处 理 方法 、 效 用 函数 和 终止 条 件 等 ， 提 高 选 星 算法 的 性 能 ， 满 足 用 户 定 位 的 精度 或 者 实时 
性 需求 。 
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图 7GPS+BDS+GLO 不 同 截止 高 度 角 的 选 星 结果 比较 
Fig. 7Comparison of GPS+BDS+GLO satellite selection results at elevation mask angles 5°, 15°and 30°. (al- 
3) number of satellites selected; (b1-3) GDOP values; (c1-3) computing time. 
4 结论 


随 着 多 全 球 卫星 导航 系统 飞速 发 展 ， 可 见 的 导航 卫星 数目 日 益 增 多 ， 接 收 机 的 性 能 日 益 增 强 ， 选 星 数 
目 也 成 为 选 星 问题 中 重要 的 考量 因素 。 本 文 将 多 系统 选 星 问题 作为 有 约束 条 件 的 多 目标 优化 问题 进行 研究 ， 
并 给 出 了 多 目标 选 星 问题 的 数学 描述 。 这 种 多 目标 选 星 方式 ， 可 以 同时 优化 几何 精度 因子 和 选 星 数目 这 两 
个 目标 ， 这 样 可 以 同时 满足 接收 机 定位 解 算 的 精度 和 实时 性 需求 ， 与 传统 的 固定 数目 选 星 方法 相 比 ， 减 小 


了 先 验 性 需求 ， 具 有 更 好 的 灵活 性 和 机 动 性 。 
在 此 基础 上 ， 本 文 提出 了 基于 NSGA-II 算法 的 多 目标 选 星 方法 。 利 用 选 星 问题 的 序 贯 性 ， 使 用 上 一 时 
刻 的 选 星 结果 生成 下 一 时 刻 的 初始 解 集 ， 大 大 减少 计算 的 复杂 度 ， 提 高 了 选 星 的 效率 ， 这 样 的 改进 方式 也 


可 以 运用 于 其 他 的 选 星 方法 中 。 同 时 改进 了 约束 处 理 方式 ， 通 过 保留 一 定 的 不 可 行 解 保持 种 群 的 多 样 性 ， 
防止 局 部 收敛 ， 提 高 了 算法 的 边缘 搜索 能 力 。 进 一 步 结合 选 星 问题 的 特征 ， 选 择 了 合适 的 遗传 算 子 和 用 来 
决策 的 效用 函数 ， 保 证 了 算法 的 可 用 性 。 通 过 仿真 证 实 了 方法 的 实时 性 和 有 效 性 ， 能 够 在 降低 接收 机 负担 
的 同时 获得 良好 的 定位 精度 ， 而 且 选 星 过 程 不 依赖 卫星 的 几何 位 置 分 布 ， 特 别 适 用 于 有 障碍 或 者 遮挡 的 复 


杂 情 况 。 
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Multi-objective and Fast Satellite Selection MethodBased on 
the NSGA-II Algorithm 
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Abstract: For multi-Global Navigation Satellite Systemconstellations, satellite selection can be considered as a 
constrained multi-objective optimization problem to optimize the Geometric Dilution of Precision (GDOP) and 
number of selected satellitessimultaneously, in order to reduce computation of the receiver and improve positioning 
accuracy. Here we describe a method based on the NSGA-II algorithm, which utilizes the sequentialcharacteristic of 
satellite selection to generate the initial population, as well as improved constraint handling, appropriate genetic 
operatorand utility function to achieve fast satellite selection.The simulation resultsshow that the proposed method 
has high reliability and real-time performance and does not rely on the geometric distribution of satellites, making 
it suitable for complex environments with obstacles and shelters. 
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